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The search for an intermediate mass Higgs boson at the LHC needs the quantitative under- 
standing of the two-photon background. A calculation of two-photon production in hadronic 
collisions at full next-to-leading order is described. It includes photons originating from the 
hadronization of QCD partons which play an important role at the LHC. A prediction for the 
invariant mass distribution for photon pairs at the LHC is presented and finally the residual 
scale dependencies are discussed. 



1 Introduction 

The understanding of electroweak symmetry breaking is a major motivation for present and 
future collider physics. Inside the Standard Model (SM) the LHC is expected to detect the 
predicted Higgs boson, if its mass is below say 800 GeV. On one hand the direct search at 
the LEP experiments place a stringent lower bound whereas on the other hand the precision 
measurements indicate an upper bound for the mass of the Higgs bosonLU: 



Also for the light Higgs boson of the minimal supersymmetric standard model a lower bound of 
88 GeV exists, where on the contrary the upper bound of around 130 GeV is dictated by the 
structure of the interaction terms. The quartic couplings are related to the gauge couplings here, 
not a free parameter as in the SM. The most promising signal for such neutral Higgs bosons 
in the mass range of 80 GeV < Mniggs < 140 GeV comes from the gluon fusion process with a 
subsequent decay of the Higgs boson into a photon pair. Though the branching fraction is small, 
B(H -> 77) - 2 x 10~ 3 in the SM, one expects a sharp peak in the invariant mass distribution 
of the photon pair, because of the narrow Higgs width of a few MeV. Though the signal is very 
sharp one is confronted with a huge background. To get a quantitative understanding of the 
signal to background ratio one has to know the background as best as one can. 
The background can be split into three contributions: 

Direct: Both photons originating from the hard partonic interaction. 

Fragmentation: At least one of the photons is created in the hadronization of a QCD parton. 
Meson decay: The photons are decay products of mesons like tt°, rj, and so on. 

Whereas the latter is essentially reducible the former are commonly called irreducible back- 
grounds though one can reduce the fragmentation contribution substantially by imposing iso- 
lation criteria, i.e. one restricts the hadronic energy inside a cone (defined in rapidity and 
azimuthal angle space) around the photon momentum to be less than some value. 



107.7 GeV < M Hiqqs < 188 GeV 




(1) 



In the following the calculation of this background at full next-to-leading order is shortly 
sketched. The inclusion of the NLO fragmentation contributions is a step beyond the calculations 
of@£lcl. A prediction of the background M 77 distribution at LHC is presented and its stability 
under variation of the different scales entering the computation is discussed. 

2 Photon pair production in hadronic collisions at NLO 

The Born term for photon pair production, qq — > 77, is of order a 2 . The radiative corrections 
consist out of all virtual and real emission corrections which lead to an order a 2 a s correction. 
One observes that in subprocesses like qg — > 77c/ a final state singularity is present if the photon 
is collinear to the quark. To absorb these singularities into a photon fragmentation function D 7 / q 
one has subsequently also to take into account other partonic reactions like e.g. qg — > jq. As 
D^j q behaves like a/a s in the investigated kinematic regime, the power counting of the couplings 
is the same as the one of the Born term. As both photons can be collinear to external partons also 
both photons should be allowed to be created in hadronization. Two-fragmentation processes 
are needed for consistency together with the respective higher order corrections, seei for more 
details. Finally there is another sizable contribution, gg — > 77, where the gluons/photons are 
attached to a fermion loop. Formally a higher order correction, this contribution is not negligible 
due to the large gluon flux at the LHC. 

The diagrammatic calculation depends on three unphysical scales: 

fi: Renormalization scale due to ultraviolet divergence, 

H-fact'- Factorization scale due to QCD initial state collinear singularities, 

M frag '• Fragmentation scale due to final state collinear singularities. 

The relation between the hadronic and the partonic cross sections for a fragmentation process 
is done by folding the partonic cross sections with the /Uj^CT-dependent parton distribution 
functions and the /i^^G-dependent photon fragmentation functions -C^/jejq,^}- 

3 The M 77 distribution at the LHC 

The comparison of our calculation with recent Tevatron data is showing a good agreement i. 
This is especially true for infrared save quantities as the M 77 distribution In Fig. [l] we plot 
the invariant mass distribution of a photon pair at the LHC with and without applying isolation 
cuts. As an example we demand the transversal hadronic energy to be less than Exmax = 5 GeV 
in a cone R = 0.4 around the photon defined in rapidity and azimuthal angle space through 
R = \J {y — y-j) 2 + {<P — <P-y) 2 - Apart from these, standard cuts were taken as indicated in the 
figure. The scale choice is = [ifact = H-frag = M 77 /2. For all figures the MRST2 set of 
parton distribution functions H and the BFG photon fragmentation functions 13 were used. To 
get an idea of the residual scheme dependence of our result for the lower curve in Fig. [l| (with 
isolation), we varied the renormalization and factorization scales and compared to that curve by 
plotting [da/dM^(n,^ FACT ) - da/dM 11 (M- yl /2, M 77 /2)]/[da/dM 77 (M 77 /2, M 77 /2)] in Fig. |. 
Due to isolation the fragmentation scale dependence is marginal 1 and is kept equal to ufact for 
simplicity. In Fig. |2| one observes that there is an accidental stability if ^,^fact are changed 
in the same directions from // = 11 fact = M 77 /2 to /jl = [if act = 2M 77 . This is due to the 

"Infrared sensitive observables are showing also good agreement for the tails of the distributions. In the infrared 
sensitive region resummation of large logarithms is needed 0. 

6 This is because we consider the fragmentation contributions in the calculation. Thus the result is of beyond 
leading logarithmic accuracy. 
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Figure 

Applying no photon isolation criterion leads to the upper curve whereas the lower curve shows the effect of 
isolation. Whereas in the upper curve the one- fragmentation contribution is dominant, in the lower curve mainly 

the direct component is present. 



fact that at LHC for the given parton energies the effect of larger (smaller) [i - means smaller 
(larger) NLO corrections - is compensated by the increased (decreased) gluon flux if [ifrag is 
also increased (decreased). By varying in an anti-diagonal way one gets a more conservative 
estimation of the uncertainty. One observes variations between 10 and 20 per cent. 



4 Conclusion 



A calculation for the photon pair production at hadron colliders at full NLO exists and is 
implemented into a partonic event generator, called DIPHOxl. Tevatron data are well described 
by the code and predictions for the LHC show a residual scale dependence of the order of 10 to 
20 per cent ||. For infrared sensitive observables, like e.g. Pt distributions, resummation has to 
be included in the calculation. In contrary the M 77 distribution is now known at NLO not only 
in the direct but also in the fragmentation part which allows more reliable quantitative studies 
of effects due to isolation criteria. 
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c To get a further stabilization at least the dominant processes would have to be calculated one order higher. 
E.g. the corrections to the box contribution have to be included but also the corrections to the gluon induced 
subprocesses. 
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Figure 2: The behaviour of the lower curve of Fig. |lj under scale variations. The relative difference is plotted for 
diagonal and anti-diagonal variations of n and (Ifact- The fragmentation scale is chosen as /if rag = [ifact- 
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